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Abstract 

We study the effect of bulk fermions on electroweak precision observables in a recently proposed 
model with warped extra dimensions and no custodial symmetry. We find that the top-quark mass, 
together with the corrections to the ZbiJji vertex and the one-loop contribution to the T parameter, 
which is finite, impose important constraints that single out a well defined region of parameter space. 

o : 

New massive vector bosons can be as light as ~ 1.5 TeV and have large couplings to the tR quark, 
and suppressed couplings to the ti-, b^ and lighter quarks. We discuss the implications for searches of 
^D \ models with warped extra dimensions at the LHC. 
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I. INTRODUCTION 
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Models with warped extra dimensions [lfl are excellent candidates for physics beyond the 
Standard Model. They solve the hierarchy problem, provide a rationale for the observed flavor 
structure, and can potentially give striking signatures at the LHC. Besides, they are weakly cou- 
pled duals of strongly coupled four-dimensional (conformal) theories and can therefore provide 
an intuition and quantitative predictions for models of strong electroweak symmetry breaking 
(EWSB). 

Models based on a purely AdSs background and a minimal field content are strongly con- 
strained by electroweak precision tests (EWPT). In particular, large contributions to the Peskin- 
Takeuchi T parameter (2| can easily push the scale of new physics beyond the reach of the 
LHC 3]. The large top mass requires the third generation SU(2)l quark doublet (as well as 
the top singlet) to be relatively close to the infrared (IR) brane. When the light generations 
are localized near the ultraviolet (UV) brane, this can lead to large corrections to the ZbiJ)L 
coupling, which also impose important restrictions. A common solution invokes a custodial 



K 



symmetry [4j, and a benchmark model of anarchic warped extra dimensions was proposed in 
Ref. J5[ (see [6| for alternatives). In that model, the S parameter -together with the one-loop 
fermion contributions to T and the ZbiJjL coupling [7|, which are calculable in such models- 
imposes the strongest constraint: the new vector bosons should be heavier than 3 — 4 TeV. 1 The 
couplings of the quarks to these Kaluza-Klein (KK) vector bosons in the previous benchmark 
model were taken to be 

9t R ~ 5gsM , 9t L ,9b L ~ 9sm , 9 q ~ ~9sm/5 , (1) 

where q represents any of the light Standard Model (SM) fermions and gsM represents a typical 
SM couplings. The conclusion of LHC reach studies is that the gauge KK resonances are difficult 
to see at the LHC, mainly due to their large width and reduced production cross section, the 
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latter resulting from the small coupling to valence q uarks U U- Lucknv, this benchmark mode, 

comes with another characteristic signature: light vector-like fermions (fermion custodians) [10], 



that are much easier to see at the LHC [ll| (see 12] for an overview). 







It has been recently suggested [13] that a departure from pure AdS 5 near the infrared brane 



1 If one gives up on the anarchy assumption, the constraints can be significantly relaxed 



can substantially decrease the T paramet er ( and, although by a smaller amount, also the S 



parameter, as was previously observed in 14], ll5|]). The authors of 13] proposed a model in 



which the reduction of the T and S parameters is so effective that new vector bosons with 
masses below 1 TeV are compatible with the EWPT without the custodial symmetry. The 
analysis in [13J assumes, however, that all fermions are localized on the UV brane. This is an 
excellent approximation for the light generations but not for the third one. 

We extend the previous study by incorporating bulk fermions and describing their effects 
on EWPT once the top and bottom masses are reproduced. This is important for a number 
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of reasons: first, without the custodial symmetry the ZhJ)^ coupling is not protected J4J and 
can impose significant constraints on the model; second, the loop-level contributions to the 
EW observables are expected to be strongly dependent on the localization of the third family 
of quarks, and can also impose significant constraints; finally, the collider implications of the 
model depend crucially on the localization of the third generation quarks. Nevertheless, we find 
that the signatures associated with the region of parameter space favored by the EW precision 
constraints are still characterized by Eq. (JTJ, although lighter resonances than in the AdSs 
background may be allowed. 

The outline of the paper is as follows. We introduce the model in Section [Til and discuss 
bulk fermions in Section IIII1 The EW constraints of the model are investigated in Section IIVI 
and its collider implications in Section [V] We conclude in Section I VI I 

II. A WARPED MODEL WITHOUT CUSTODIAL SYMMETRY 

The model under consideration represents a departure from AdSs in the infrared. We simply 



summarize in this section the main results 



'or the gravitational background and the bosonic 



field content, while referring the reader to 13] for full details. We will then introduce bulk 



fermions in the model and discuss their effects on the EWPT. The gravitational background is 
given by the metric 



with warp factor 



ds 2 = e-' 2A{y) r 1 ^dx ll dx 1 ' - dy 2 , (2) 



A(y) = ky-\\og(l-^) (31 



The extra dimension is bounded by two branes: the UV brane, localized at y = 0, and the IR 
brane, localized at y — y\. The warp factor would vanish at y s but we will always consider 
parameters such that this metric singularity remains hidden behind the IR brane (y s > yi). 

The gravitational parameters are therefore u, j/i, y s and k (the curvature scale at the UV 
brane, of the order of the Planck mass) which is assumed to set the scale of all dimensionful 
5D quantities. We will trade y s for the value of the curvature radius at the IR brane, given in 
units of k by 

kin = ^^ ~ Vl) . (4) 

^1 _ 2z/2/5 + 2 ^k(y s - yi ) + v^k\y s - y x f 



Requiring perturbativity of the gravitational expansion bounds its value by kl\ > 0.2 13|. 
The position of the IR brane, y\, can be fixed by requiring the gravitational background to 
generate the Mp/TeV hierarchy. We will simply set A(yi) = 35, which determines y\. 

The bosonic content of the model consists of the SM gauge fields, with Neumann boundary 
conditions (BC) on both branes, and an electroweak scalar doublet, the Higgs. Other scalars 



involved in the gravitational background and the stabilization of the interbrane distance 13J are 
irrelevant for this discussion, although their phenomenology, in particular that of the radion, 
can be interesting. The bulk gauge bosons can be expanded in KK modes (we focus on the \i 
component here, which is the relevant one for EWPT and collider implications): 

M*,v) = ^Y,fn{y)^\x), (5) 



/j/i 

where the profiles satisfy 

[d y e- 2A d y + ml] f* = , d y f A \ y=0m = , (6) 

(we treat EWSB perturbatively), and are normalized according to 

- / dyf A fi = 6 nm . (7) 

The boundary conditions fix the value of the KK masses. They are of the order of the effective 
IR scale 

k cS = A'( Vl ) e~ A ^ , (8) 

i.e. of order the warped down curvature at the IR brane. We show in Fig. [1] the mass of the 
first gauge KK mode, in units of k e g, for different values of v and kl\. 
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FIG. 1: Mass of the first gauge KK mode in units of the effective IR scale k e s, defined in Eq. ([8]), as 
a function of v and for different values of kL\, 



For the Higgs field we can write 



H(x,y) = -^-e ix{x > y) [ 

V2 \h(y)+Z(x,y) 



where the Higgs KK expansion reads 



(9) 



(n). 



X 



(10) 



and the wavefunctions are normalized to 



1 f yi 
- / dye- 2A f n f m = 6 nm . 

yi Jo 



(ii) 



For a light Higgs we can assume that the Higgs vacuum expectation value (vev) is carried by 
the zero-mode, which has a profile 



/«(y) « N h e ak y, 



with Nh fixed by the normalization condition ( II ip . Following [13] we trade a for 



(12) 



-2(a-2)ky s 



ky s [-2(a-2)ky s 






T[l--,-2(a-2)k(y s - yi [ 



(13) 



which is a measure of how much fine-tuning in the 5D parameters we need to impose in order 
to preserve the Randall- Sundrum solution to the hierarchy problem [1]. 



In summary, we have two free parameters for the gravitational background, v and hL\\ and 
one for the Higgs background, 5. However the EWPT are quite insensitive to the latter so we 
will simply fix 5 = 0.1 in the following. 

III. BULK FERMIONS 

We introduce now bulk fermions in the model. The quadratic part of the action for bulk 
fermions reads 



S 



fd A x dy e~ AA tf \e A i0 + (d 5 - 2A')-f 5 - M~\ V , (14) 



where M is a 5D Dirac mass allowed by the symmetries. Following the standard procedure we 
expand the bulk fermions in KK modes 



Vi 

where we have defined the 4D chiralities ^ l,r — |(1 : F7 5 )$, and the KK profiles satisfy the 
following system of equations: 

[d y - (2A> + M)] J* = -e A m n ti , 

[d y - (2A! - M)\ ti = e A m n f« , (16) 

with the ortho-normality conditions given by 

i rvi i rvi 

- / dye- 3A titi = ~ / dye- 3A f*f r R n = 5 nm . (17) 

Vi Jo 2/i Jo 

Bulk fermions with a left-handed (LH) zero-mode are obtained by imposing Dirichlet boundary 

conditions for their right-handed (RH) chirality on both branes. Such bulk fermions will be 

denoted as "[++] fermions". The boundary conditions for the opposite, LH, chirality are fixed 

by the equations of motion. Similarly, bulk fermions with a RH zero-mode are obtained by 

imposing Dirichlet boundary conditions on their LH chirality, and will be denoted as "[ ] 

fermions". The LH (RH) massless profile of a [++] ([ ]) bulk fermion is given by 

tf>*(y) = n^^SWW), (18) 
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FIG. 2: Value of the localization parameter, C1/2, that makes a LH fermion zero-mode mostly delo- 
calized, as a function of v and for different values of kL\. 

where N ' is fixed by the normalization condition. In the following we will consider a constant 
bulk mass Mi = Qk, where % denotes the fermion type. 2 

In the present metric background there are no exactly flat fermion solutions, unlike for a pure 
AdSs background, ^4rs(|/) = ky, with a constant 5D Dirac mass defined by c = 1/2. However, 
the c parameter does control the localization of the fermion zero-modes, either towards the IR 
or UV brane. In particular, there is always a background-dependent value, ci/2, such that for 
c < Ci/2, the (LH) fermion zero-mode is mostly IR localized, while for c > c\i<i it is mostly UV 
localized. This value plays a role analogous to c = 1/2 in the pure AdS 5 background. In Fig. [2] 



we show the value of Ci/ 2 , defined such that e 2j4 '°Vo'(0) 



e 2 



Mvi) fL 



ft(y x ) (see Eq. (El and 



footnote |5]) , as a function of the input parameters v and hL\ (see also Fig. |7J). 

Let us now discuss the couplings of the fermion zero-modes to the Higgs and gauge boson 
KK modes in this model. The gauge couplings are given by 



£ D 95 dye~ SA tyA^ 



V— = dye 
^ Jyi J 



LfL fA 



fLfL £ 
3A J n J mJ r 

Vi 



i> 



(n) 



A i (m) 



(L^R) 



T,9^?* r ¥r ) +(L^R) 



(19) 



2 In the limit in which y s — > y\ there could be problems of strong coupling similar to the ones present in 
soft- wall models [15|, and a y— dependent mass term might be necessary [16j (see [17| for other realizations 
of flavor in soft- wall models). 



whereas the Yukawa couplings can be computed from 
C D Y 5 fdy e~ AA Q L H U R + (L^ R) (20) 

v , f Q(n) , (0) ,U(m) 



where we have defined /i = ^ < -°- ) (x). The 5D gauge coupling can be fixed by matching the 
coupling of the gauge zero-mode to the observed 4D coupling. Assuming a tree-level matching 
this gives 

95 = \/yi9i ■ (21) 

For the 5D Yukawa coupling, Ref. [18] finds, based on NDA [l9|, that its maximum value is 
given by 

Y 5 < Y 5 max « 4^ » ( 22 ) 



which corresponds to strong coupling at a scale of the third KK level. We plot in the left panel 

of Fig. [3]the Yukawa coupling between the zero-modes of a [++] fermion Q, and a [ ] fermion 

T, assuming Y5 = Y 5 max , as a function of a common localization parameter c = cq = —ct- In the 
right panel we show the coupling of a LH fermion zero-mode to the first two gauge boson KK 
modes divided by the coupling to the gauge boson zero-mode (assuming tree level matching) 
as a function of c. We have chosen v = 0.4, and kL\ = 0.2 in both plots. We see the change 
of behavior around C1/2 ~ 0.93, which corresponds to the most delocalized zero-mode fermion 
profile (see Fig. EJ) . 

For c > C1/2, we see that the Yukawa coupling becomes exponentially suppressed. Further- 
more, in the same region the coupling to the gauge boson KK modes becomes almost universal 
(i.e. independent of c). Thus, as for the pure AdSs case, the assumption that the light fermions 
are exactly localized on the UV brane is an excellent approximation (as far as EWPT are 
concerned) in these models. 

IV. ELECTROWEAK CONSTRAINTS 

If flavor is explained by means of localization in the extra dimension, the effects of new 
physics in models with warped extra dimensions are almost universal and can therefore be en- 



coded in the oblique parameters of 20(, with the most important exception being the couplings 
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FIG. 3: Left panel: Yukawa coupling, A = Aqq , as a function of the bulk mass for cl, = — cr = c and a 
5D Yukawa coupling that saturates the maximum value given in Eq. (|22p . Right panel: coupling of a 
LH fermion zero- mode to the first two gauge boson KK modes in units of go = g^l \[V\- I n both cases 
we have taken v = 0.4 and kL\ = 0.2, which lead to maximally delocalized fermions for c ~ 0.93. 

of the bottom quark. This can be most easily seen with holographic methods. However, the 
collider implications of these models are easier to understand if we discuss how the corrections 
are generated in the physical basis (i.e. the KK-basis). A general discussion of EWPT in mod- 
els with warped extra dimensions and the equivalence of different methods to compute them 



is given in 



2l| . Here we follow the equivalent notation in [13| to make the comparison easier. 



Before presenting the detailed EW precision analysis, we comment on the calculability of such 
effects in the present class of models. 

A. Calculability 



As already emphasized, the present class of models does not have a custodial symmetry to 
protect the Peskin-Takeuchi T parameter, nor the corrections to certain gauge-fermion cou- 
plings. As a result, it is possible to write a term in the bulk that violates the custodial 
symmetry, £ 5 D (k/A 3 ) \H^DmH\ 2 , where we wrote the coefficient in units of the 5D cutoff, 
A. The dimensionless factor, k, is UV sensitive. However, when the Higgs propagates in the 
bulk, the 1-loop contributions are finite by power-counting. To see this, it is simplest to use a 
normalization where the gauge bosons, Wm, have mass dimension 1 [we write the kinetic term 
as (— l/4pf) Wf rIN W aMN }. Then, the contributing diagrams are simply proportional to g\ or 
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FIG. 4: Examples of radiative corrections to the oblique parameters, as well as to the non-oblique 
vertex corrections. The most important contributions arise from the top KK-tower in the loop. Upper 
row, left to right: one-loop contributions to T, S and 5gb L , respectively. Lower row: examples of the 
corresponding two-loop contributions. 

yf, where both the 5D (top) Yukawa coupling, y 5 , and the 5D gauge coupling have mass di- 
mension —1/2. (We show in the upper-left corner of Fig. H]an example diagram with a fermion 
loop.) It follows that the piece proportional to 77^ in the loop integral has mass dimension — 1, 
and therefore it is IR dominated. Subtracting the zero-mode contribution, the remainder is 
controlled by the scale k e g. At two-loop order, an example of which is shown in the lower-left 
corner of Fig. HJ the diagram is logarithmically divergent by power counting. 

A similar point has been made in connection to the 5-parameter when the Higgs is taken 
as a bulk field. 3 The middle diagram in the upper row of Fig. H] is finite, while at 2-loop 
order a logarithmic UV sensitivity is encountered (e.g. in the middle diagram of the lower row 
of Fig. Hj), corresponding to the bulk operator H^a a HW^, IN B MN . Similar remarks apply to 
"vertex corrections" of the form (HW M H)($T Ml $) or (H^a a D M H)(^T M a a ^). Examples of 
1- and 2-loop contributions to these vertex corrections are shown in the right column of Fig. HI 
We also note that operators localized on the branes (corresponding either to T, S or 5gb L ) can 
be induced only at three (and higher-order) loop order. 

We conclude that just allowing the Higgs to propagate in the bulk can make the oblique 
parameters effectively calculable: the incalculable pieces associated with the 5D local operators 
above are suppressed compared to the finite, one-loop contribution roughly by M KK /A. In the 



K. Agashe, private communication. See also 
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particular models studied here, the localization properties of the Higgs, to be discussed in the 
next subsection, imply that the couplings of the Higgs to the KK fermion or gauge states are 
suppressed compared to the situation in an AdS 5 background. In fact, for the favored region of 
parameter space such couplings are well in the perturbative regime (thus resembling more closely 
the case of flat rather than typical warped extra-dimensions). However, although higher-order 
contributions due to Yukawa or weak gauge couplings can be expected to be further suppressed, 
there remains some uncertainty associated with higher order QCD contributions. As we will see 
in detail below, one of the effects of the departure from AdSs of the gravitational background 
is to push the KK modes closer to the IR brane. This effect increases the coupling among KK 
modes, thus reducing the scale of strong coupling in the QCD sector. This is reminiscent of 
the position dependent cut-off in soft- wall models [14j, [l5| , and affects the EW observables at 
two-loop and higher order. 4 Nevertheless, we find that the finite one-loop contributions to the 
EW observables -especially to T- can be significant, and therefore one should not neglect such 
radiative effects. To get a concrete idea about their impact on the EW fit and the resulting 
bounds on the KK scale, we will assume that the QCD strong coupling scale is high enough to 
make the 5D description a reasonable approximation. Furthermore we will assume that 2-loop 
and higher-order QCD effects do not dramatically change the one-loop results. One should, 
however, keep in mind that QCD effects may not be negligible. 

B. The Oblique Corrections at Tree Level 

We perform now a detailed analysis of the EW precision constraints in the class of models 
without custodial symmetries under discussion. In this subsection we focus on the oblique 



analysis at tree level (which was already performed in 13|]), and in the next one we take the 
most important non-oblique contribution into account, i.e. the correction to the Zbi}>L coupling, 
as well as the one-loop contributions. This will allow us to better understand the impact of 
the various effects. Let us start with the gauge KK-modes. The effects of new physics can 
be classified in three types: corrections to the gauge boson self-energies, to fermion-gauge 



4 A smaller hierarchy in the spirit of the Little RS model [22] would reduce the tension with the low QCD 
cut-off and increase the coupling of the KK-gluon to light quarks, thus having an important impact on collider 
phenomenology. 
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(23) 



FIG. 5: Contribution to the coefficients a, (3 and 7 of Eq. (|24p . The double line represents the tower 
of massive gauge boson KK modes, the red squares the mixing between the gauge KK and zero-modes, 
and the blue dot the coupling of the light fermions to the gauge KK modes. 

couplings and to four-fermion interactions. Each of these effects is characterized by coefficients 
denoted by ft, (3 and 7, respectively. For the present case, only gauge bosons obeying Neumann 
boundary conditions on both branes are relevant, in which case one has 13] 

fy ye ^(n h {y)-^\ , 

[ Vl dye^(n h (y)-y-)({l f (y)-y 
Jo V Vi/ \ Vi 

The function Q(y) is denned by 5 

o m - l /V f ~\ n - / e ' 2A{y) ^(y) 

tti{y) = — / dyujiiy), Ui{y) = < 

where fi t o(y) is the wave function for the scalar or fermion zero- modes (see previous sections). 
With our normalization we have Qi(yi) = 1. These coefficients are diagrammatically shown in 
Fig. [5l The wavy double line represents the tower of massive gauge boson KK modes, resummed 
in Eqs. ( 1241) . The blue dot represents the coupling of the SM fermions to the new vector bosons, 
and the red square the mixing between SM and heavy gauge bosons through Higgs insertions. 
The oblique parameters due to the KK physics can be written in terms of these coefficients as 



scalars 
fermions 



(24) 



T 

S 

W 



g' 2 v 2 



(ft-2/3 + 7) 



2 

sV(-/3 + 7 ) 



Y 



g 2 v 2 



"7 



(25) 
(26) 
(27) 



5 The u>i(y) are nothing but the (square of the) "physical wavefunction" profiles, i.e. the profiles with warp 
factors taken out, as dictated by dimensional analysis (which just redshift all mass scales appropriately). 
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FIG. 6: 95% CL lower bound on the mass of the first KK gauge boson, obtained from a tree-level 
oblique analysis. We assume UV localized fermions and different values of the input parameters v and 
kL\. See text for details on the fit procedure. 

where g and g' are the SU{2)l and U{1)y couplings, respectively, and v = 174 GeV is the 
Higgs vev. These are related to the Peskin-Takeuchi S and T parameters [2| by aS = As^yS 
and aT = T. We also include a Higgs contribution to S and T, given by [23] 

1 



AS 
AT 



2tt 



[g{m 2 h /m 2 z ) - g(m 2 Tei /m 2 z )] 



[f (ml/ml) - fimlJml)] 



16vrc^ 



where 



g(y) 



f(y) 



dxx(5x — 3) ln(l — x + yx) 



In cly-lny 



+ 



Iny 



(28) 
(29) 

(30) 
(31) 



cw - y c w( l - y) 

and we take m re f = 115 GeV in the fit. 

We show in Fig. [H]the result of a fit to the oblique parameters when all fermions are assumed 
to live on the UV brane for different values of the input parameters v and kL\. We have used 
an updated version of the code in 24J, obtaining the bounds as follows. For fixed values of v 
and kLi (i.e. for a fixed metric background), we compute the minimum of the x 2 ? including 
only the Z-pole observables, 6 as a function of the Higgs mass (imposing the direct LEP bound 



6 These are the 26 observables associated with the W mass (two measurements), the Z-line shape and lep- 



13 



>2 4 



CS 




0.0 



0.2 



1,11,11111,11,11,11 



0.4 



0.6 



y/ji 



0.8 



1.0 



FIG. 7: Fermion zero-mode, Higgs and first gauge KK mode profiles, y/uj(y) for v = 0.4 and kL\ = 0.2. 



of ~ 114 GeV) and the value of the IR scale k e g. We then fix the Higgs mass to its value at the 
minimum and compute the bound on k e g by requiring A% 2 = 3.84 (95% CL for one degree of 
freedom). The resulting value of the first gauge boson KK mode mass (see Fig. [I]) is plotted as 
a function of v for different values of kL\. We see that masses as light as ~ 1 TeV are allowed 
for kL\ = 0.2 and v ~ 0.45. This plot reproduces the results presented in 13] up to small 
differences (< 100 GeV), due to the different fit procedure. 

The reason for the reduced lower bound on the mass of the gauge KK modes, compared to 
the pure AdSs background, despite the absence of custodial symmetry in the model, is easy 
to understand. In the AdSs case, the warp factor forces the gauge boson KK modes to be 
localized close to the IR brane, where the Higgs is also localized. The large overlap then makes 
the mixing between the gauge boson zero-mode and massive KK-modes large (the red square 
in Fig. [5]). The light fermions on the other hand are localized towards the UV brane and 
therefore their coupling to the gauge boson KK modes (the blue dot) is typically small. This 
is the reason for the usual enhancement of the T parameter (a is proportional to the gauge 
mixing squared) over the S parameter, and of S over W and Y . The departure from the AdSs 
background has several effects that go in the right direction to improve the EWPT. The first 



ton FB asymmetries (8), heavy flavor (6), effective sin 2 dw (2), and 8 leptonic/strange quark polarization 
asymmetries. We do not fit the SM parameters, but only those associated with the new KK physics. 
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FIG. 8: Couplings of UV localized fermions to the first-level KK gauge bosons, in units of the zero- 
mode gauge coupling. We show curves as a function of v, and for different values of kL±. 

is that the KK gauge bosons are more strongly localized towards the IR brane whereas the 
Higgs, although localized towards the IR brane, reaches a maximum before the IR brane [the 



"physical" Higgs wavefunction is given by ^ojh(y) ~ e~ A ^ e aky , and the suppression near the 
IR brane arises from the nearby singularity in A(y)]. This reduces, sometimes dramatically, 
the gauge mixing through the Higgs vev, and therefore the contribution to the T parameter 
(and somewhat less the 5* parameter). This effect can be observed in Fig. in which we show 
the "physical" profiles y/ui of Eq. (|24l) . for the Higgs, the first gauge KK mode, and for bulk 
fermion zero-modes with three different values of the bulk masses. The maximum of the Higgs 
profile before the IR brane is evident, together with the very strong localization of the first 
gauge KK mode, leading to a reduced value of the relevant overlap integral. 

The second effect of the warp factor is that the coupling of UV localized fermions to gauge 
KK modes is reduced with respect to the standard AdS 5 case. This further reduces the S 
parameter (but will also have a negative impact on collider searches). In Fig. [8] we show the 
value of the coupling of UV localized fermions to the first gauge boson KK mode (in units of 
9o — 95/ 'x/yT) f° r different values of v and kL\. Only for small values of v and large values of 
kLi does one approach the coupling of the AdSs background. Note that for the values preferred 
by the EWPT (see Fig. [6]) this coupling is reduced to almost half the AdS 5 value. 
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C. The Effects of the Third Generation 

Let us now consider the effect of the third quark generation. We consider three bulk fermion 
fields with the following quantum numbers under the SU(2) L x U(l)y gauge group and BC 

Q=(2,l/6)~[++], T = (1,2/3) ~[— ], B = (l,-l/3)~[— ] , (32) 

with localization parameters cq, ct and cb, respectively. For fixed values of these localization 
parameters we can choose the 5D Yukawa couplings so that the top and bottom masses are 
reproduced. Due to the upper bound on Y 5 , the top mass cannot be generated if the Q and T 
zero- modes are far from the IR brane. The LH bottom, which is in the zero-mode of Q, can 
then receive large corrections to its couplings, from gauge and fermion KK modes. Anomalous 
contributions to the ZbiJ)L coupling can therefore impose stringent constraints in warped models 
without custodial protection. The general expression for the tree-level correction to the ZbiJjL 



coupling induced by the gauge boson KK modes is given in 2l|. For our model it reads 



g 2 v 2 



2 6 



0Q ~ Puv) , (33) 



where the term proportional to fiuy corresponds to the universal part that is absorbed in the 
oblique parameters. To this we have to add the tree-level fermionic contribution, that we have 
computed exactly by diagonalizing numerically the KK fermion mass matrix, and computing 
the resulting coupling to the Z of the lightest mass eigenstate (the bottom quark). We will 
consider the 1-loop effects separately. 

In order to test the dependence of the constraints on the assumptions on the 5D Yukawa 
coupling we have taken the following benchmark scenarios. For each value of v and kL\ and 
each value of cq we fix ct so that the top mass is reproduced (including the effect of the mixing 
with fermion KK modes) assuming that Y£ = y^ max (scenario 1). We note that, due to the 
maximum of the Higgs profile, there is a fixed value of ct for which the overlap is maximal 
and therefore the 5D Yukawa coupling is minimal. Therefore, we also consider the case that 
Yg is the minimal one for which it is possible to reproduce the correct top mass (scenario 2). 
Regarding the bottom sector, we consider three different scenarios by fixing the 5D bottom 
Yukawa to Y" 5 6 = Y$ (scenario a), Y$ = Y 5 4 /5 (scenario 6), Y" 5 6 = Kjf/10 (scenario c). 
Scenario a assumes exact anarchy whereas in the other two we allow for deviations between 
different Yukawas. 
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FIG. 9: Tree-level 95% CL lower bound on the mass of the first gauge KK mode as a function of 
the localization of the LH top/bottom multiplet (cq) for the six scenarios discussed in the text. We 
fix v = 0.4 and kL\ = 0.2. The lines end at the value of cq beyond which the top mass cannot be 
generated. For comparison we also show the result of the fit if the effects of third generation quarks 
are neglected (horizontal "ST" line). The maximal derealization is obtained for c ~ 0.93 (see Fig. [2]). 



We have studied the tree-level effect of third generation quarks on EWPT by performing a 
scan over v and kL\. For each value of these parameters, we have computed the maximum value 
of cq that allows to generate the top mass (i.e. the furthest from the IR brane). We have then 
scanned over the values of cq smaller than this maximal value for the six different scenarios 
(la, 16, lc, 2a, 26, 2c) described above. The fit has been performed as in the case of the gauge 
contribution, described in the previous subsection, but including now the constraint from the 
Zb[J)L coupling described above. Fixing cq and marginalizing over the Higgs mass, we find the 
value of k e s that gives A% 2 = 3.84, corresponding to 95% CL for one degree of freedom. In all 
cases the preferred value of the Higgs mass is close to its current lower limit rrih ps 114 GeV. 
We show in Fig. [9] a sample result for v = 0.4 and kL\ = 0.2. We display, based on a tree-level 
analysis, the 95% CL lower bound on the mass of the first gauge KK mode as a function of cq 
for all six scenarios, and also the result of the fit when all fermions are localized on the UV 
brane as considered in [13| (horizontal line denoted by ST in the plot). It is clear that the bound 
is very sensitive to the value of the localization of the LH top/bottom doublet as expected for 
a non-custodial model. The correction is smaller for a larger cq (the LH top/bottom further 
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FIG. 10: One-loop contribution to T, S and 6gb L , as a function of cq, in scenario 2c (see text). We 
take v = 0.4, kL\ = 0.2 and k c ^ = 1 TeV (which corresponds to Mkk ~ 2.3 TeV). 

from the IR brane). The non-trivial result we find is that the top mass (just barely) allows cq 
to be large enough to make the corrections to the Zbibi coupling negligible (the lines in the 
figure end at the point of cq beyond which the top mass cannot be generated). Unfortunately, 
loop effects change this conclusion, as discussed below. The other property that is clear from 
the plot is that exact anarchy (scenarios la and 2a) is extremely constrained by EWPT. The 
reason is that the Yukawa couplings in the bottom sector are very large and the mixing with 
the bottom KK modes induces very large corrections to the ZolOl coupling. However, we see 
that a suppression in the 5D bottom Yukawa by a factor of 5 is already enough to get the bound 
on the scale of new physics reasonably low. We have also checked that the i^^bjiW^ coupling 
is in these cases typically < 10~ 4 (in units of g/y/2) and should therefore cause no trouble with 
b — > sj constraints 25] (a full analysis, including the fermion KK modes in the calculation, to 
properly account for this constraint, is beyond the scope of this paper). 

As mentioned in Subsection IIV A\ the 1-loop contributions to T, S and 5gf, L are finite (while 
at two- loop order, they are logarithmically divergent). We have computed these one-loop effects 
using the methods described in [7|, and show a numerical example in Fig. [lOj Here we have 
taken a gravitational background with v = 0.4 and kL\ = 0.2, and we have assumed scenario 2c 
defined above, with k e g = 1 TeV. Recall that scenarios 2 are such that the 5D Yukawa coupling 
is the minimal one that still allows to reproduce the top quark mass, for given cq. In this sense, 
these scenarios minimize the size of these one-loop effects, which are controlled by this coupling. 
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FIG. 11: Left panel: 95% CL lower bound on the mass of the first gauge KK mode as a function of 
the localization of the LH top/bottom multiplet (cq) in scenario 2c, including one-loop effects. The 
different lines correspond to different values of rrih, with the one marked as "m/j at Xmin" corresponding 
to marginalization over m/j. The curves are terminated where the goodness-of-fit gives a 5% likelihood. 
Right panel: comparison of the tree-level oblique, full tree-level, and tree-level plus 1-loop bounds on 
Mkk, assuming marginalization over m^. In all cases, we fix v = 0.4 and kL\ = 0.2. 

We see in the figure that the 1-loop contribution to the T parameter can be significant, and 
imposes an upper bound on cq that is stronger than the one coming from the top mass itself. By 
contrast, the 1-loop contributions to S and 5gt, L are relatively small (in the figure we show 20 x S 
and 10 3 x Sgb L )- This constraint is in tension with the one due to the (tree-level) modification 
of the Zhibi coupling. We show in the left panel of Fig. [11] the resulting lower bound on the 
KK-gluon mass corresponding to scenario 2c. The thick solid (blue) line corresponds to the 
fit procedure used at tree-level, i.e. evaluating A% 2 = 3.84 with a Higgs mass that minimizes 
the total x 2 (typically near the LEP Higgs bound). However, we also show the bounds on 
the KK scale assuming other fixed values of the Higgs mass (as would be appropriate if a 
Higgs of such a mass was actually discovered). We see from the thick solid blue curve that 
-marginalizing over m^- a lower 95% CL bound of Mkk ~ 2.3 TeV is found for cq ~ 0.88. We 
note that for this c Q one has Xmm/ dof = 25.6/24 at (m h , M KK ) = (114 GeV, 4.3 TeV), which 
gives a goodness-of-fit with 37% likelihood. This results mainly from a compromise between 
the (tree-level) 5gt, L and the 1-loop contribution to T. However, the black dot-dashed line 
shows that for a heavier Higgs, with fixed m/, = 500 GeV, a KK-gluon as light as 1.5 TeV 
(1.4 TeV) would be possible, which corresponds to the 95% CL contour about a best fit point 
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with Xmin/dof = 34/25 (x 2 min /doi = 37.6/25) leading to a 10% (5%) likelihood. Thus, it may 
be possible to have warped models with KK-gluons around 1.5 TeV that fit the EW data 
reasonably well. As a summary plot that highlights the impact of the various contributions 
to the EW observables discussed above, we show in the right panel of Fig. [TT] the bounds on 
the KK-gluon mass, as a function of cq, for the EWPT fits at a) tree-level in the oblique 
approximation [thin solid black line], b) full tree-level [dotted brown line] and c) tree level plus 
1-loop [thick solid blue line]. In all cases we marginalize over rrih, although as just pointed out 
this may lead to an overly pessimistic conclusion in regards to how low Mkk could actually be. 
At any rate, it is clear from this figure that both the corrections to the Zbibi coupling and the 
one- loop effects play an important role in determining the allowed Mkk- 

We also point out that the models consistent with the EWPT up to 1-loop order, which have 
a relatively low Mkk, always have KK-fermion Yukawa couplings that are perturbative. For 
instance, at the minimum of the solid thick blue curve in Fig. [Til with Mkk ~ 2.3 TeV, the 4D 
Yukawa couplings of the form hQ£ Ir are all 0(1), while the off-diagonal ones (i.e. coupling 
different KK levels) are much less than one (becoming smaller the further apart the masses of 
the two KK modes). 7 This is a result of the suppression in overlap integrals associated with the 
non-trivial profile of the Higgs field, much as in the KK-gauge/Higgs couplings illustrated in 
Fig. [7J Thus, higher order (divergent) effects involving additional powers of Yukawa couplings 
are expected to be suppressed. The most important effects that remain are associated with 
QCD higher-order corrections, as mentioned in Subsection II V Al Thus, the above results should 
be taken as an illustration of how light the KK resonances can reasonably be, as far as the EW 
precision constraints are concerned. 

These results have important implications for collider searches. First, the absence of custo- 
dial symmetry implies a quite minimal spectrum of massive modes. In particular, very light 
fermions 



10 



111 ], natural in custodial models are not expected in the model under study. We 
show in Fig. [12] the mass of the first fermion KK mode of a [++] field as a function of the 
localization parameter c for different values of v and kL\. These figures, together with Fig. [TJ 



7 By contrast, scenarios 1 have diagonal 4D Yukawa couplings for the KK fermions of order 3 — 4, and always 
lead to a very large 1-loop contribution to T, unless Mkk is above O(10 TeV). In such cases, higher-order 
contributions associated with the Yukawa interactions may not be suppressed, and can have a large impact 
on the EW analysis. Nevertheless, barring tuned cancellation between these and the UV contributions, one 
expects that the KK resonances will be out of the LHC reach in such scenarios. 
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FIG. 12: Mass of the first fermion KK mode for a [++] field in units of k e s as a function of its 
localization parameter (with respect to the mostly delocalized value) for values of v = 0.3, 0.4, 0.5, 
0.6, 0.7 (top to bottom, left panel) and kL\ = 0.2, 0.3, 0.4, 0.5, 0.6 (bottom to top, right panel). 

show that, in some cases, the first KK gluon can decay into a top (or bottom) and a heavy 
fermion, although the latter are in this case much closer to threshold than in custodial models. 
Channels with a heavy fermion of the first and second generations are also typically open: un- 
der the anarchic assumption such decay channels are very suppressed, but they could be more 
important in other scenarios for flavor. Second, masses lighter than previously considered in 
models with warped extra dimensions (with semi-anarchic Yukawas) may be allowed by EWPT. 
This result is very sensitive to the localization of the third generation quarks which determine 
the coupling of the top and bottom quarks to the gauge KK modes, as well as the KK-fermion 
Yukawa couplings. We show in Fig. [13] the couplings of the different quarks to the first KK 
gluon in units of go = g^J ■s/Vi for scenarios 2, with v = 0.4 and hL\ = 0.2. We see that the 
largest is the coupling to t R i R , followed by QlQl and then to light quark pairs. Thus, the KK- 
gluon decays dominantly into RH top pairs (or, perhaps, a channel involving one KK-fermion 
and the associated zero-mode), bearing some resemblance to the scenario of Eq. (JTJ. However, 
the reduced couplings to the light quarks can make the discovery more challenging, for a given 
M KK - 
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FIG. 13: Couplings to the first gauge KK mode (in units of go = g^j \fy\) of the third generation 
LH quarks Ql, RH top £#, RH bottom 6r, and light quarks q, as a function of cq, for scenario 2 
(sub-scenarios a, b and c give similar results for the couplings). We fix v = 0.4 and kL\ = 0.2. The 
localization parameters ct, cb and the top and bottom masses have been fixed as described in the 
text. The bn and q couplings are almost identical. 

V. PHENOMENOLOGICAL IMPLICATIONS 

In this section we present the main collider implications of some selected points in parameter 
space. As we have seen in the previous section, the EW fit prefers models of type 2, which 
always have the KK-gluons relatively strongly coupled to t R (see Figs. [11] and [13]) . Thus, the 
ti decay channel for the KK-gluon is always dominant. To illustrate the expected signal at the 
LHC, we consider two models, as suggested by the analysis in the previous section. 

The first one corresponds to scenario 2c with v = 0.4, kL\ = 0.2 and cq = 0.88. With 
rrih = 114 GeV, we find a x 2 /dof = 25.6/24, giving a likelihood of 37%. This fit is only slightly 
worse than the SM one: the total x 2 is slightly reduced, but there is one additional degree of 
freedom, corresponding to k e g. The KK gluon mass is 



M 



KK 



2.3 TeV 



(95% CL) 



(34) 



We also find that the first KK resonance of the third generation quark SU(2) L doublet has a 
mass Mq ps 2.1 TeV, while the first KK resonance of the bottom SU(2)l singlet has a mass 
M B w 1.85 TeV. Both are sufficiently light for the decays G^ -> QfCffl and G^ -> b { §B$ 
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to be open. The corresponding couplings to fermion pairs are 

g Qh w 1.95& , g tR « 4.63# s , g bR w 0.02& , # g w -0.13& , (35) 

V.Q^ w 3 - 85 &» ' 9 bR:B m « 1-1% S , (36) 

where we omitted the superscripts for the zero modes, and the first line refers to SM fermion 
pairs. The first KK-resonance of the top SU(2)l singlet is heavier than the KK-gluon, so that 
this channel is kinematically closed. However, the first KK resonances of all the remaining 
SM fermions are lighter than the KK-gluon, and are therefore open as decays of the form 
qW _). g(°)g( 1 ). Nevertheless, the relevant couplings are all significantly smaller than g s , so that 
these are somewhat rare decays that, in spite of the multiplicity, do not change appreciably the 
width of the KK-gluon. With the above couplings, we find that r G (i) ~ 710 GeV, and that the 
KK-gluon has the following branching fractions: 

BR(G^ -► ft) « 0.81 , BR(G^ -»■ 66) w 0.12 , BR(G^ -» qq) « 0.004 , (37) 
BR(G^ ->• t { l ] t) « 0.02 , 5i?(G« ->■ &£ ) 6) « 0.04 , 5i?(G« -)■ 6^6) « 0.01 , (38) 

where in the last line both conjugate processes (e.g. G^ 1 ' — > t L t and G^ 1 ) — >■ t^ t) are under- 
stood. The total G^ production cross section is ~ 24 fb (~ 188 fb) at a CM energy of 7 TeV 
(14 TeV). 

The second model has instead cq = 0.906 which, with m^ = 500 GeV, leads to x 2 /dof = 
34/25, giving a likelihood of 10% (still reasonably large). The KK-gluon mass is now 

M KK « 1.5 TeV (95% CL) , (39) 

while the first KK resonance of the third generation quark SU(2)l doublet has a mass Mq ps 
1.3 TeV, and the first KK resonance of the bottom SU(2)l singlet has a mass Mb ~ 1.2 TeV. 
The corresponding couplings to fermion pairs are 

g QL » 1.17 , g tR « 4.62 , g bR » -0.03& , g q » -0.13& , (40) 

Now we have r G <i) ps 390 GeV and the following branching fractions: 

BR(G^ ->■ ft) = 0.83 , BR(G^ ->■ 66) = 0.05 , BR(G^ ->■ gg) = 0.005 , (42) 

BR{G^ ->• 4 1} t) « 0.03 , £i?(G« -► 6^6) « 0.06 , £#(£« -► 6^6) « 0.02 . (43) 
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FIG. 14: ti invariant mass distribution in the SM (red dotted), in the model with extra-dimensional 
physics (solid blue) and the contribution of just the KK gluon exchange (dashed purple). The left 
panel corresponds to a KK-gluon with Mkk ~ 1-5 TeV and the couplings of Eqs. (|40p - (|4ip . The right 
panel corresponds to Mkk ~ 2.3 TeV and the couplings of Eqs. ([35]) - ([36]) . 



The total G*- 1 -* production cross section is ~ 0.19 pb (~ 1.3 pb) at a CM energy of 7 TeV 
(14 TeV). 

n 

The BR's in these models are relatively similar to the benchmark of Ref. [5|, but with 
lighter masses and reduced couplings to the light quarks. Note, however, that the first model 
has a non-negligible branching fraction into bottom pairs, and that there are "exotic" channels 
involving a KK fermion with BRs at the few percent level (in the case of the third generation; 
for the first two generations the corresponding BRs are expected to be much smaller, although 
the precise values depend on the details of how flavor is implemented). 



We have implemented these models in MADGRAPH/MADEVENT v4 [26J, using PYTHIA 



2811 for detector simulation. We use the 



6 [23] for hadronization and showering and PGS4 
CTEQ6L1 parton distribution function, with the QCD renormalization and factorization scales 
equal to the central m? of the event. We have then implemented the ti CMS analysis of Ref. [29J , 
which we use to normalize our SM sample, generated by MG/ME, once we check that the shape 
of the ti invariant mass distribution is reasonably reproduced. We use this (re)normalization 
factor for the signal + background, taking into account interference effects. We consider here 
the semi-leptonic ti channel, separating the e and \x channels 29]. We show in Fig. [14] the 
expected ti invariant mass distribution at the LHC, at the partonic level. The left (right) panel 
corresponds to the second (first) model above at \/s = 7 TeV (y/s = 14 TeV). In both plots we 
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represent the SM prediction with a red dotted line and the prediction of the model (including 
the interference with the SM it contribution) in solid blue. Also, just to guide the eye, we show 
the contribution assuming only the KK gluon exchange as a purple dashed line. Although the 
lightest mass case shows a slight excess over background, these results suggest that extracting 
the signal will be challenging. 

VI. CONCLUSIONS 



We have studied a recently proposed model in five dimensions with a warped background 13] 
that can be consistent -from the point of view of EW precision constraints- with new gauge 
boson resonances at about 1.5 TeV. This represents a significant reduction of the bound on 
the new physics scale compared to fairly elaborate models based on AdSs backgrounds. The 
scenario is also relatively simple, with just the SM field content promoted to 5D, plus a new 
"stabilizing" scalar field that is responsible for the all-important deviation from AdSs near 
the IR brane. A simple way to understand the relaxation of the bound can be obtained by 
considering the localization of the relevant fields: the gauge boson KK modes are strongly 
attracted towards the IR brane, while it is still possible to accommodate a profile for the Higgs 
field that, although mostly IR localized (so that the RS solution to the hierarchy problem is 
preserved), attains its maximum before reaching the IR brane. As a result, any overlap integral 
that involves the Higgs and the (first) gauge or fermion KK modes can be dramatically reduced 
compared to cases where all the fields are localized within a distance of order 1/k from the IR 
brane, as happens in the AdSs background. Since the deviations from the SM are determined 
by such overlap integrals (and the scale of the new resonances), lighter states can be allowed. 

Here we have extended the study in Ref. [13] to the case where the SM fermions can be arbi- 
trarily localized along the extra-dimension, which allows to understand the SM flavor structure 
as arising from the localization of fermion fields. We point out that reproducing the observed 
top quark mass, which requires the top quark to be localized sufficiently close to the IR brane, 
can result in significant constraints from anomalous contributions to the ZbiJjL coupling, be- 



yond the constraints from the oblique parameters considered in Ref. 13|. This is a direct 
consequence of the absence of any custodial protection j4j] together with the "anarchical" ex- 
planation of flavor. In addition, we find that the 1-loop corrections to the EW observables, 
which are finite in these scenarios, strongly restrict the allowed region of parameter space. In 
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particular, a tension between the (tree-level) corrections to the ZbJ)L coupling and the 1-loop 
contribution to the Peskin-Takeuchi T-parameter, strongly constrains the localization of the 
third generation quarks and the resulting 95% CL bound on the KK scale. Nevertheless we find 
that for an optimal top/bottom SU(2)l doublet localization, and when the 5D bottom Yukawa 
coupling is a factor of a few smaller than the top one, the KK-gluons can be as light as 1.5 TeV 
if the Higgs is heavy, with the likelihood inferred from the EW fit still being reasonably large. 
However, if the Higgs mass is allowed to float in the fit, one finds a 95% CL lower bound of 
M KK ~ 2.3 TeV, 

A KK-gluon mass of about 1.5 TeV can open the exciting possibility that the Randall- 
Sundrum solution to the hierarchy problem, based on warped compactifications, and with a 
bulk SM field content, could lead to observable resonances at the LHC with lower luminosities 
than previously thought. Compared to the widely studied AdSs framework, the result arises 
from two opposing effects: the previously mentioned strong localization of the gauge boson 
resonances towards the IR brane implies a reduction of its couplings to the light fermions, thus 
leading to a suppression in production. This can be compensated by the lower allowed mass of 
the gauge KK modes. However, we find that discovering such a resonance in the dominant ti 
channel is likely to be challenging. Boosted top techniques and a very detailed knowledge of the 
ti tail will likely be required to discover these modes. We also find that, unlike in the case of an 
AdSs background without the custodial symmetry, the present class of models allows for the 
production of single KK-fermion resonances of the SM fields, in KK-gluon decays. Although 
these decay channels are likely subdominant (in particular, they do not dramatically change 
the KK-gluon width), they may provide an interesting handle on the fermionic resonances. 8 
Nevertheless, rather high luminosities are likely necessary. One should recall, however, that we 
are working under the (semi-) anarchic assumption of flavor. If the light families are closer to 
the IR brane (with a correspondingly smaller 5D Yukawa coupling) the production cross section 
can be larger, and relatively light spin-1 resonances in the 1-2 TeV range could be more readily 
observed at the LHC. 



Note that in models with custodial symmetry, fermion custodians are expected to be in general lighter and 



these channels are likely to be more relevant [32|, |33j. The discovery of new fermions with exotic charges, like 



Q = 5/3 1CJ, |ll|, would in any case be a clear smoking gun of custodial models. 



2(3 



Acknowledgments 

We would like to thank Kaustubh Agashe, Mert Aybat, Cedric Delaunay, Gero von Gersdorff 
and especially Manuel Perez- Victoria for useful discussions. A.C. and J.S. are supported by 
MICINN (FPA2006-05294,FPA2010-17915, FPU and Ramon y Cajal programs) and Junta de 
Andalucia (FQM 101, FQM 03048 and FQM-6552). E.P is supported by DOE grant DE-FG02- 
92ER40699. 



[1] L. Randall, R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370-3373. hep-ph/9905221]; Phys. Rev. 



Lett. 83 (1999) 4690-4693. |hep-th/9906064|. 
[2] M. E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381 (1992). 

[3] S. J. Huber, C. -A. Lee, Q. Shah, Phys. Lett. B531 (2002) 112-118. | hep-ph/01 11465] . 
[4] K. Agashe, A. Delgado, M. J. May, R. Sundrum, JHEP 0308 (2003) 050. | hep-ph/0308 036 1 ; 

K. Agashe, R. Contino, L. Da Rold, A. Pomarol, Phys. Lett. B641 (2006) 62-66. 

[hep-ph/0605341| . 
[5] K. Agashe, A. Belyaev, T. Krupovnickas, G. Perez and J. Virzi, Phys. Rev. D 77, 015003 



(2008) [arXiv:hep- ph/0612015|; B. Lillie, L. Randall and L. T. Wang, JHEP 0709 (2007) 074 
|arXiv:hep-ph/0701166| ; A. Djouadi, G. Moreau, R. K. Singh, Nucl. Phys. B797, 1-26 (2008). 
[arXiv:0706.4l9Tl [hep-ph]]. 
[6] M. S. Carena, A. Delgado, E. Ponton, T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 68 



(2003) 035010 arXiv:hep-ph/0305188|: Phys. Rev. D 71 (2005) 015010 |arXiv:hep-ph/0410344|; 



A. Djouadi, G. Moreau and F. Richard, Nucl. Phys. B 773 (2007) 43 [arXiv:hep-ph/0610173] . 



[7] M. S. Carena, E. Ponton, J. Santiago and C. E. M. Wagner, Nucl. Phys. B 759 (2006) 202 
[arXiv:hep-ph/0607106] ; Phys. Rev. D 76 (2007) 035006 [arXiv:hep-ph/0701055| . 

[8] C. Delaunay, O. Gedalia, S. J. Lee, G Perez, E. Ponton, Phys. Rev. D83, 115003 (2011). 
[arXiv: 1007.02431 [hep-ph]]: la7Xiv:1101.2902l [hep-ph]. 

[9] A. L. Fitzpatrick, J. Kaplan, L. Randall, L. -T. Wang, JHEP 0709 (2007) 013. | hep-ph/070il~50~l ; 



K. Agashe, H. Davoudiasl, G Perez, A. Soni, Phys. Rev. D76 (2007) 036006. | jhep-ph/0701186| ; 
K. Agashe et al, Phys. Rev. D 76 (2007) 115015 [arXiv: 0709 .00071 [hep-ph]]; K. Agashe, 
S. Gopalakrishna, T. Han, G Y. Huang and A. Soni, Phys. Rev. D 80 (2009) 075007 



27 



[arXiv:0810.l497l [hep-ph]]; K. Agashe, A. Azatov, T. Han, Y. Li, Z. G. Si and L. Zhu, Phys. 
Rev. D 81 (2010) 096002 |arXiv:091 1.00591 [hep-ph]]. 

[10] G. Cacciapaglia, C. Csaki, G. Marandella and J. Terning, Phys. Rev. D 75 (2007) 015003 
[arXiv:hep-ph/0607146] ; R. Contino, L. Da Rold and A. Pomarol, Phys. Rev. D 75 (2007) 
055014 [arXiv:hep-ph/061204 8 1 ; A. D. Medina, N. R. Shah and C. E. M. Wagner, Phys. Rev. 
D 76 (2007) 095010 [arXiv: 0706 .12811 [hep-ph]]: A. Pomarol and J. Serra, Phys. Rev. D 78 (2008) 
074026 |arXiv:0806.3247l [hep-ph]]. 

[11] R. Contino and G. Servant, JHEP 0806 (2008) 026 |arXiv:0801.1679l [hep-ph]]; A. Atre, 
M. Carena, T. Han and J. Santiago, Phys. Rev. D 79 (2009) 054018 |arXiv:0806.3 966 [hep- 
ph]]. J. A. Aguilar-Saavedra, JHEP 0911 (2009) 030 [ arXiv:0907.3155l [hep-ph]]; J. Mrazek and 
A. Wulzer, Phys. Rev. D 81 (2010) 075006 |arXiv:0909.3977 [hep-ph]]; F. del Aguila, A. Car- 
mona and J. Santiago, JHEP 1008 (2010) 127 [arX iv: 1001.5151 [hep-ph]]; Phys. Lett. B695 
(2011) 449-453. [arXivjl00L|206] [hep-ph]]; G. Dissertori, E. Furlan, F. Moortgat, P. Nef, JHEP 
1009 (2010) 019. |arXiv:1005.44l4l [hep-ph]]: A. Atre, G. Azuelos, M. Carena, T. Han, E. Ozcan, 
J. Santiago, G. Unel, [arXiv: 1102. 19871 [hep-ph]]. 

[12] J. Santiago, PoS ICHEP2010 (2010) 557. [ arXiv: 1103.4114 [hep-ph]]. 

[13] J. A. Cabrer, G. von Gersdorff, M. Quiros, Phys. Lett. B697 (2011) 208-214 [ arXiv:1011.2205l 
[hep-ph]]; JHEP 1105 (2011) 083 [a rXiv: 1103. 1388 I I [hep-ph]]: |arXiv:1104.3l"49l [hep-ph]]. 

[14] A. Falkowski and M. Perez- Victoria, JHEP 0812 (2008) 107 |arXiv:0806. 17371 [hep-ph]]: JHEP 



0912 (2009) 061 [ arXiv:0901.3 777 [hep-ph]]. 
[15] B. Batell, T. Gherghetta and D. Sword, Phys. Rev. D 78 (2008) 116011 [arXiv: 0808 .39771 [hep- 

ph]]; 



[16] S. Mert Aybat and J. Santiago, Phys. Rev. D 80 (2009) 035005 [arXiv:0905.3032 [hep-ph]]; AIP 



Conf. Proc. 1200 (2010) 611 [arXiv:0909.3999 [hep-ph]]. 
[17] T. Gherghetta and D. Sword, Phys. Rev. D 80 (2009) 065015 [arXE:0907.3523 [hep-ph]]; A. Del- 



gado and D. Diego, Phys. Rev. D 80 (2009) 024030 [arXiv:09 05.1095 [hep-ph]]; M. Atkins and 
S. J. Huber, Phys. Rev. D 82 (2010) 056007 |arXiv:1002.5044l [hep-ph]]. 



[18] K. Agashe, A. Azatov, L. Zhu, Phys. Rev. D79, 056006 (2009). [arXiv:0810.1016 [hep-ph]] 



[19] Z. Chacko, M. A. Luty and E. Ponton, JHEP 0007, 036 (2000) ]arXiv:hep-ph/9909248|. 

[20] R. Barbieri, A. Pomarol, R. Rattazzi and A. Strumia, Nucl. Phys. B 703 (2004) 127 



28 



[arXiv:hep-ph/0405040] . 
[21] H. Davoudiasl, S. Gopalakrishna, E. Ponton, J. Santiago, New J. Phys. 12 (2010) 075011. 

[arXiv: 0908 .19681 [hep-ph]]. 
[22] H. Davoudiasl, G. Perez, A. Soni, Phys. Lett. B665 (2008) 67-71. [a rXiv: 0802 .02031 [hep-ph]]. 
[23] M. J. G. Veltman, Acta Phys. Polon. B8 (1977) 475. 
[24] Z. Han, W. Skiba, Phys. Rev. D71 (2005) 075009. |hep-ph/0412166] ; Z. Han, Phys. Rev. D73, 

015005 (2006). |hep-ph/0510l25] . 



[25] P. L. Cho, M. Misiak, Phys. Rev. D49, 5894-5903 (1994). |hep-ph/9310332|; K. Fujikawa, A. Ya- 



mada, Phys. Rev. D49, 5890-5893 (1994); G. Burdman, M. C. Gonzalez-Garcia, S. F. Novaes, 



Phys. Rev. D61, 114016 (2000). |hep-ph/9906329|. 
[26] J. Alwall et al., JHEP 0709 (2007) 028 |arXiv:0706.2334l [hep-ph]]. 
[27] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605 (2006) 026 [arXiv:hep-ph/0603175| . 



[28] PGS4 http://www.physics.ucdavis.edu/~ conway/research/software/pgs/pgs4-general.htm 



[29] CMS collaboration Note: TOP-10-007-PAS. 

[30] A. Abdesselam, E. B. Kuutmann, U. Bitenc, G. Brooijmans, J. Butterworth, P. Bruckman de 

Renstrom, D. Buarque Franzosi, R. Buckingham et ai, EPHJA,C71,1661. 2011 C71 (2011) 

1661. |arXiv:1012.54T2l [hep-ph]]. 
[31] B. Lillie, J. Shu, T. M. P. Tait, Phys. Rev. D76, 115016 (2007). [ arXiv:0706.3"960l [hep-ph]]. 
[32] M. Carena, A. D. Medina, B. Panes, N. R. Shah and C. E. M. Wagner, Phys. Rev. D 77 (2008) 

076003 |arXiv:0712.0095l [hep-ph]]. 



N. Vignaroli, arXiv: 1107.4558 [hep-ph]. 



29 



